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多 变量 追踪 研究 的 模型 整合 与 拓展 : 
考察 往复 式 影响 与 增长 趋势 


х] 源 
(西南 大 学 心理 学 部 , 认 知 与 人 格 教育 部 重点 实验 室 , 重庆 400715) 


摘 要 追踪 研究 当中 , 交 又 滞后 模型 可 以 探究 多 交 量 之 间 往 复式 影响 ， 潜 增长 模型 可 以 探究 个 体 增 长 趋 
势 。 对 两 类 模型 进行 整合 , 例如 同时 关注 往复 式 影响 与 个 体 增长 趋势 AN TARLM SRA, MAES Ж 
异 成 分 ,衍生 出 随机 蕉 距 交 叉 滞 后 模型 、 特 质 - 状 态 -误差 模型 、 自 回归 潜 增 长 模型 、 结 构 化 残 差 潜 增 长 模型 
等 。 以 交 又 滞后 模型 和 潜 增 长 模型 分 别 作 为 基础 模型 ,从 个 体 间 /个 体内 交 蜡 分 解 的 角度 对 上 述 各 类 模型 梳理 ， 
整合 出 此 类 模型 的 分 析 框 架 ， 并 拓展 建立 “因子 结构 化 潜 增 长 模型 (factor latent curve model with structured 
reciprocals)” 作 为 统合 框架 。 通 过 实证 研究 (早期 儿童 的 追踪 研究 -幼儿 园 版 , ECLS-K), 建立 21049 名 儿童 的 阅 
读 和 数学 能 力 的 往复 式 影响 与 增长 趋势 。 研 究 发 现 , 分 离 了 稳定 特质 的 模型 拟 合 最 优 。 研究 也 对 模型 建 模 思 
路 和 模型 选择 提供 了 建议 。 

关键 词 ”追踪 研究 ,往复 式 影响 ,增长 趋势 ， 因 子 结构 化 潜 增 长 模型 
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追踪 研究 (ongitudinal analysis) 也 叫 纵向 研究 、 描述 总 体 增长 趋势 。 对 于 均值 差异 的 比较 ,主要 
纵 贯 面 研究 , 是 在 比较 长 的 时 间 内 , 对 一 个 个 体 或 ”考虑 到 消除 系统 偏差 、 减 小 误差 方差 等 因素 ,可 
一 些 个 体 的 心理 发 展 进行 有 系统 的 定期 的 研究 , 在 以 用 重复 测量 的 方差 分 析 或 多 元 方差 分 析 等 方 
心理 学 、 教 育 学 、 社 会 学 、 医 学 等 学 科 中 被 广泛 使 法 。 第 二 类 研究 旨 在 探讨 多 变量 之 间 的 往复 式 影 
用 。 随 着 当代 统计 技术 的 不 断 发 展 ,特别 是 社会 科 响 (reciprocal relations, 或 相互 影响 ),， 与 更 严格 的 
学 统计 方法 的 发 展 , 运用 结构 方程 模型 、 多 层 线 性 因果 影响 (causal effecb 进 行 区 分 ( 温 忠 麟 , 2017)。 
模型 等 统计 方法 对 追踪 数据 进行 研究 ,能 关注 更 多 可 以 通过 定义 交叉 滞后 模型 (Cross-Lagged Model, 
И == Бл 


的 信息 ， 比 如 增长 趋势 的 个 体 差异 ， 以 及 影响 个 CIM), 探讨 多 变量 变化 的 相互 影响 。 最 后 一 类 研 


А» 


体 差异 的 各 层级 水 平 的 因素 等 。 由 此 可 见 ， ә: 究 关注 的 是 更 为 宏观 的 增长 趋势 或 发 展 轨迹 ， 引 


研究 越 来 越 细 化 ， 所 关注 的 问题 也 越 来 越 深入 。 | 进 漠 增 长 模型 (Latent Growth Model, LGM). 
结合 以 往 的 研究 (Usami et al., 2019; 刘 红 云 ， 在 长 期 的 应 用 当中 ， 对 上 述 几 个 问题 的 研究 


ЈЕ, 2003; FEE, “ЖНЖ, 2019), 可 以 总 结 出 往往 都 是 孤立 的 。 随 着 近期 一 些 模型 的 产 出 ， 研 


追踪 研究 主要 关注 三 个 方面 的 问题 : (1) 均 值 差 异 究 者 可 以 按照 研究 问题 ， 同 时 或 选择 性 地 考察 往 
的 比较 ，(2) 对 多 变量 相互 影响 进行 解释 ,以 及 (3) 复式 影响 和 增长 趋势 。 从 应 用 着 的 角度 ， 整 合 此 
类 的 追踪 研究 将 要 解决 的 问题 ,给 研究 者 提出 新 


的 研究 框架 和 研究 建议 。 
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后 介绍 结合 了 两 个 基础 模型 的 整合 模型 。 
11 交叉 滞后 模型 : 多 变量 往复 式 影 响 

探讨 变量 之 间 的 往复 式 影响 常用 的 模型 是 交 
又 灌 后 模型 (Cross-Lagged Model, CLM), 又 称 为 
果 模 型 、 交 叉 清 后 面板 模型 (Cross-Lagged Panel 
Model)、 自 回归 交叉 滞后 模型 (Autoregressive Cross- 
Lagged Mode]) 等 (Kenny & Harackiewicz, 1979); 后 
来 随 着 结构 方程 模型 (Structural Equation Modeling, 


变化 和 外 在 偶然 因素 对 其 变化 的 影响 。 对 于 传统 
CLM 模型 一 般 满 足 2 次 的 测量 就 能 达到 识别 要 求 
(Usami et al., 2019); 而 广义 交叉 滞后 模型 则 需要 至 
少 4 次 测量 才能 识别 (Zyphur Allison, et al., 2020)。 
研究 者 对 CLM 模型 中 最 感 兴 趣 的 模型 参数 
为 自 回归 系数 8 和 交叉 滞后 回归 系数 ( 刘 文 等 ， 
2015)。 此 类 模型 揭示 的 变量 之 间 的 关系 通常 被 研 
究 者 定义 为 “往复 式 影 响 (reciprocal effects)”, БІ 


SEM) 技 术 的 发 展 ， 能 够 很 容易 地 被 研究 者 所 应 用 
(例如 J6reskog，1970)。 交 叉 滞 后 模型 以 自 回归 为 
基础 (公式 1), 分 析 所 测 变量 的 跨 时 稳定 性 
(temporal stability, 参见 Hamaker et al., 2015)。 


Xit = Myt 十 Xy (1) 
其 中 ， Ж. = Ж.” жа, о A 回归 模型 中 ， 观 
测 指标 x 表示 被 试 i 在 第 t 次 测量 上 的 值 ， 系数 pt 


随 着 时 间 的 推移 ， 每 个 变量 的 变化 都 会 影响 其 他 
变量 的 变化 (Usami et al., 2019)。 但 值得 注意 的 是 ， 
往复 式 影响 和 因果 影响 仍然 不 能 划 等 号 ， 后 者 需 
要 满足 更 多 的 前 提 假 设 。 一 般 地 , 在 严格 的 实验 
设计 下 ， 如 果 已 经 控制 了 无 关 变 量 ，CLM 可 以 用 
来 揭示 变量 之 间 的 因果 关系 ( 温 忠 记 , 2017)。 

сім 模型 定义 相对 单 ， 对 测量 次 数 要 求 也 不 多 ， 
也 能 在 一 定 程度 上 揭示 不 同 变 量 随时 间 变 化 的 相互 


为 自 回归 系数 ， 描 述 了 变量 本 身 随时 间 变 化 的 效 
应 大 小 ， 即 跨 时 稳定 性 ,数值 越 大 则 稳定 性 越 高 ; 
dxit 为 残 差 , па 为 方程 的 截 距 。 容 易 拓展 到 多 变量 
情形 , 例如 含有 两 个 变量 的 情况 (公式 2, 图 1)。 


ж 
Xin = Mat + Xit 


* (2) 
Vit = Hye + Vit 

其 中 ， Ж, = i) + yi жа. у = Ву + 
дуа, ТАНЯ В 和 By DA, уа 
Жуу 26 ИЕЛ AR, УХА x (或 变量 y) 在 第 
1 次 的 测量 受到 第 1-1 次 变量 y (或 变量 x) 的 影响 ; 
da 和 dy 表示 残 差 ， 服从 联合 正 态 分 布 。 有 研究 
者 (Zyphur Allison, et al., 2020) 认 为 交叉 沸 后 影响 
不 应 仅 存在 于 变量 本 身 ， 环 境 ( 变 量 之 外 的 偶然 因 
素 ) 也 可 能 对 个 体 的 后 续 表现 产生 影响 。 故 自 回归 参 
数 与 交叉 滞后 参数 拓展 到 第 六 1 次 残 差 对 第 + 次 测量 
的 影响 ,被 称 之 为 广义 交叉 滞后 模型 (General Cross- 
Lagged Model; 30) 等 , 2021)。 区 别 于 变量 的 跨 时 
稳定 性 ， 残 差 对 下 一 次 测量 的 影响 包含 了 “偶然 因 
素 ” 的 贡献 ， 使 得 模型 能 同时 处 理 变量 本 身 的 前 后 


dn 
图 1 交叉 滞后 模型 示意 图 


影响 关系 ,应 用 领域 涉及 到 个 体 特质 、 教 育成 果 与 
精神 疾病 特征 等 领域 的 长 期 性 、 滞 后 性 影响 的 追踪 
研究 (Cacioppo et al., 2017; Chiu & Du, 2019; Liu & 
Hou, 2018; Selzler et al., 2019)。 国 内 这 方面 研究 主要 
集中 在 青少年 问题 行为 、 抑 郁 、 能 力 发 展 的 影响 因 
素 等 方面 ( 陈 红 君 等 , 2019; ЖО 等 , 2018)。 
12 жака: 个 体 增 长 趋势 

另 一 类 问题 是 描述 个 体 增长 趋势 ,包含 总 体 
增长 形态 的 特点 及 其 是 否 存 在 个 体 差异 ( 即 发 展 
群体 的 同 质 性 )。 潜 增长 模型 (Latent Growth Model, 
LGM, Bollen & Curran, 2006) 是 建立 在 SEM 视角 
下 ， 通 过 定义 带 有 均值 结构 的 验证 性 因素 分 析 模 
型 ， 固 定 特殊 载荷 大 小 来 实现 增长 趋势 的 描述 。 


以 线性 LGM 为 例 ( 公 式 3~ 公 式 4， 见 图 2): 
测量 部 分 yi, = Mi +E- Di; + En (3) 
结构 部 分 - 2 a 
hi =K +6) 


其 中 , yi 表示 被 试 i 在 第 1 次 的 观测 指标 , 其 
测量 误差 分 别 对 应 st。 定 义 两 个 潜 变 量 : 截 距 Wi 
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图 2 潜 增 长 模型 示意 图 


和 线性 斜率 12i， 潜 变量 因子 的 均值 分 别 对 应 кү 和 
кә, 方差 对 应 于 var (0 和 var (02), Са 和 ;服从 联 
RESIM HoF, WEE mi 和 poi 的 系数 矩阵 A 
可 以 根据 时 间 和 形态 灵活 定义 或 增 减 斜率 因子 ， 
诸如 二 次 增长 模型 、 不 定义 增长 曲线 模型 和 单 因 
素 增长 曲线 模型 等 ( 唐 文 清 等 , 2014; Жаа, E 
向 阳 , 2018)。 除 了 连续 的 增长 之 外 , 还 可 以 定义 间 
断 增长 的 趋势 ， 如 多 阶段 增长 模型 (piecewise 
growth model, Liu et al., 2018; Lock et al., 2018)。 
一 般 对 于 两 个 因子 的 线性 增长 模型 , 3 次 测量 是 
证 模型 识别 的 最 低 要 求 ， 而 如 果 需 要 增加 增长 
子 (如 二 次 、 多 阶段 因子 ), 需要 增加 测量 次 数 。 
描述 个 体 增长 趋势 有 助 于 研究 者 把 握 个 体 发 
展 的 形态 ， 从 而 根据 不 同 阶 段 的 发 展 特征 制定 不 
同 的 教育 政策 或 干预 手段 。 如 果 需 要 从 政策 或 根 
源 上 来 干预 增长 , 模型 中 一 般 会 考虑 非 时 变 
(time-invariant) 协 变量 ,他 们 不 会 随 着 时 间 的 变化 
而 变化 ， 其 影响 也 会 较为 持久 、 滞 后。 如 学 生 的 
性 别 、 社 会 经 济 地 位 、 历 史 事件 等 相对 稳定 的 因 
素 (Diallo et al., 2017; Liu et al., 2020)。 如 果 干 预 是 
即刻 的 ， 则 时 变 (time-variant) 协 变量 更 为 精准 。 他 
们 属于 状态 类 的 变量 ,其 影响 是 及 时 的 、 变 化 的 。 
此 时 ,多 水 平 模型 (Multilevel Model, MLM) AY 
时 变 协 变量 定义 方式 也 可 以 纳入 LGM 的 研究 当 
中 (Liu et al., 2016; McCoach & Kaniskan, 2010). 


HX 


LGM 模型 定义 也 相对 简单 ， 研 究 者 如 果 有 3 
次 及 以 上 的 追踪 则 可 使 用 该 模型 建 模 。LGM 将 变 
化 分 解 为 初始 状态 和 增长 趋势 ， 可 以 宏观 地 从 两 
个 维度 来 考察 总 体 增 长 趋势 , E SEM 框架 下 定义 
也 相当 灵活 。 国 内 的 应 用 研究 近年 逐日 增多 , R 
中 在 儿童 青少年 社会 能 力 、 孤 独 感 、 自 我 概念 、 
积极 适应 的 发 展 和 老年 人 的 认 知 功能 发 展 等 方面 
( 侯 桂 云 等 , 2018; 李 彩 娜 等 , 2017; “ЫМ 等 ， 
2017; 刘 俊 升 等 , 2013; 张 晓 , 2011)。 其 中 , 也 包 
含 时 间 跨 度 较 短 的 状态 类 变量 , 如 人 际 冒 犯 、 抑 
郁 等 ( 程 刚 等 , 2016; 杨 满 云 等 , 2017)。 
13 测量 误差 与 随机 截 距 

以 上 回顾 了 两 个 “基础 模型 ”， 可 以 在 他 们 基 
础 之 上 添加 一 些 特 殊 议题 ， 将 变异 进一步 分 解 。 
以 交叉 滞后 模型 为 载体 ， 借 用 SEM 的 思路 ,考察 
每 个 指标 的 测量 误差 ,可 以 建立 “因子 交叉 滞后 模 
型 (Factor Cross-Lagged Model, FCLM)” (Jéreskog & 
Sörbom, 1974; McArdle, 2009; Usami et al., 2015, 


见 图 3). 
测量 部 分 : 
Xin = Е + Exit (5) 
Yi = Fyi + Eyit 
结构 部 分 : 


Fa = А, + BF get) + Hak tet жа,, 
Py = y+ BF aay + ЕВ жа, 

НҢ, Fea 和 Ру 表示 潜 变 量 ， 是 分 别 排除 了 
测量 误差 sit 和 гу, 之 后 的 真 分 数 。 在 ЕСІМ 模型 
中 ,交叉 滞后 系数 都 建立 在 潜 变 量 上 ， 当 测量 误 
差 为 0 时, ЕСІМ 可 以 简化 成 标准 CLM 模型 。 与 
此 同时 , ЕСІМ 当中 每 个 测量 指标 到 因子 之 间 的 
路 径 系 数 限定 为 1, 测量 误差 服从 联合 正 态 分 布 ， 


(6) 


图 3 ”因子 交叉 滞后 模型 示意 图 
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需要 至 少 3 次 测量 才能 满足 模型 识别 的 要 求 。 如 
果 每 次 的 测量 模型 为 多 指标 ， 可 以 按照 SEM 模型 
识别 的 要 求 建立 测量 模型 。 


类 模型 进行 拓展 的 重要 出 发 点 (Bainter & Howard, 
2016; Curran & Bauer, 2011; Hamaker et al., 2015; 
Voelkle et al., 2014; ЖҰП 等 , 2021), Heh, 个 体 


男 一 方面 ， 自 回归 系数 虽然 描述 了 跨 时 稳定 
性 , 但 是 经 过 足够 长 的 时 间 ， 其 原始 的 秩序 也 会 
丢失 (如 图 1 中 xx, 的 自 回 归 间 接 效 应 为 乘积 
Bx2BxsBx4， 其 值 小 于 1D)。 随 机 截 距 交 叉 滞 后 模型 
(Random Intercept Cross Lagged Model, RI-CLM) 
的 建 模 方 式 是 将 交叉 滞后 模型 增加 一 个 随机 截 距 
项 ， 用 来 描述 真正 的 非 时 变 的 特质 (time-invariant， 
trait-like stability, 参见 Hamaker et al., 2015; 
Mulder & Hamaker, 2020， 见 图 4)。 具 体 建 模 方式 
如 下 


Xit = Hyt + hai + is (7) 
Vit = Ну + Nyi + Vit 
其 中 ， A = Baay HVV + adi, у =B Yunt 
Valia- +9 os 在 这 个 模型 中 , mix; 和 yi 表示 非 时 
变 稳定 因子 ,用 以 呈现 跟 特 质 相 关 的 因素 ,固定 
其 对 于 所 有 测量 指标 的 影响 都 为 1。 这里， 区 别 于 
LGM ЖЯ, уу, 和 ту SEN 0， 其 方差 协 方 差 矩 
Е Ф 和 截 距 项 jt 和 wy 自由 估计 。 一 般 要 求 至 少 
3 次 测量 才能 满足 模型 识别 的 要 求 。 
由 于 引入 了 特质 变量 ,可 以 将 个 体 发 展 的 变 
异 分 解 为 个 体 间 变异 (between-person difference) 
和 个 体内 变异 (within-person difference), 这 是 此 


间 变 异 来 源 为 提取 所 有 个 体 在 数 次 测量 之 间 的 公 
ЖЕ, 表示 个 体 跨 时 间 的 稳定 特征 ; 个 体内 变 
异 是 每 个 个 体 在 数 次 测量 之 间 的 滞留 影响 表示 
个 体 本 次 状态 受到 先前 状态 的 影响 (Hamaker 
et al., 2015)。 在 RI-CLM 中 , 加 入 随机 截 距 的 目的 
在 于 将 稳定 的 、 不 随时 间 变 化 的 变异 真正 分 解 为 
非 时 变 的 特质 (个 体 间 变异 ) 和 跨 时 稳定 性 (个 体内 
变异 ) 因 素 。 此 时 ， 自 回归 系数 pt 和 py 被 称 为 个 
体内 滞留 参数 (within-person carry-over, Hamaker 
et al., 2015)， 亦 被 称 为 个 体 惰 性 (inertia, Mulder & 
Hamaker, 2020)。 如 果 将 RI-CLM 的 myxi 和 myi 限 
定 成 0 时 ， 就 成 为 一 个 标准 CLM 模型 ,表示 没有 
跨 时 间 的 稳定 特质 ， 即 不 存在 个 体 间 差异 。 

在 SEM 框架 下 , 很 容易 将 测量 误差 与 随机 截 
距 模 型 结合 ， 得 到 “特质 -状态 -误差 ”模型 
(Trait-State-Error, TSE, Kenny & Zautra, 1995, 
2001)。 并 且 ， 这 个 考察 单个 变量 变化 的 模型 很 快 
拓展 到 多 元 情形 (Luhmann et al., 2011; Zautra 
et al., 1995， 见 图 5)。 首 先 , TSE 模型 考虑 测量 误 
差 ， 即 每 个 题目 都 有 各 自 的 测量 信和 度 ( 公 式 5)。 分 
解 出 测量 误差 si 和 sm 后, 真 分 数 Fit 和 Қа 的 变 
异 继续 被 分 解 成 两 个 部 分 : 一 部 分 是 稳定 特质 


图 4 随机 截 距 交叉 滞后 模型 示意 图 


JF! 
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图 5 特质 -状态 


(та 和 my 的 影响 ,一 部 分 是 状态 (sit 和 F* yin) 
之 间 的 关系 (包括 自 回归 系数 Bur Al By BOSH 
系数 yx 与 yy 两 部 分 ): 

Е = Myr + Mh + BaF oat) ер YaF ыт %а, (8) 

Fr = Hy + Th yi + ВЕ) + FF sees t+ dyin 

一 般 地 ， 特 质 与 状态 的 影响 都 固定 为 1。 如 果 
假定 时 变 的 参数 ， 比 如 自 回归 参数 1、 交叉 滞后 参 
数 人 、 残 差 和 协 方差 矩阵 都 随 着 上 的 不 同 而 自由 佑 
计 ， 要 求 至 少 4 次 测量 才能 满足 模型 识别 的 要 求 。 
如 果 假 定 非 时 变 的 参数 ， 即 上 述 参 数 都 限定 为 相 
等 , 3 次 测量 即 可 满足 模型 识别 要 求 。 

TSE 中 ， 随 机 截 距 的 ma 和 mw 代表 非 时 变 因 
子 ， 即 个 体 的 特质 (trait)， 是 “个 体 间 差异 ” A 
归 和 交叉 滞后 部 分 表示 当前 的 测量 受到 前 一 次 测 
量 的 影响 , 是 “个 体内 差异 ” 此 二 者 与 RILCLM 相 
同 , 与 此 同时 ,考察 每 个 测量 指标 的 真 分 数 ， 被 称 
之 为 个 体 的 状态 (state), 不 能 被 解释 的 部 分 是 误 
差 (error)， 此 与 FCLM 相同 。 不 难看 出 , TSE 同时 
包含 了 随机 截 距 和 测量 误差 .是 将 RI-CLM 和 
FCLM 进行 了 整合 。 也 有 研究 者 将 误差 固定 为 0 
后 将 其 看 成 是 考虑 了 测量 信和 度 的 RLCLM 模型 
(Mulder & Hamaker, 2020)。 

TSE 在 后 续 发 展 中 改名 为 “稳定 特质 自 回 归 
特质 与 状态 模型 (Stable Trait Autoregressive Trait 
and State Model, STARTS)” (Kenny & Zautra, 2001; 


误差 模型 示意 图 


Schmitt & Steyer, 1993)。STARTS 在 命名 上 区 分 了 
两 种 “特质 *: 稳定 特质 指 来 自 于 随机 截 距 的 部 分 
(ты 和 wy)， 即 “个 体 间 变 异 ” 自 回归 特质 (8 和 
px0 即 “个 体内 滞留 影响 ”。 状 态 指 每 次 测量 的 真 分 
(Frit 和 yt)， 即 排除 “测量 误差 ”所 以 在 命名 上 ， 
STARTS 更 为 精确 。 由 于 STARTS 和 TSE 只 是 命 
名 上 的 区 分 , 为 了 简化 , 仍然 将 此 模型 称 之 为 “ 特 
质 - 状 态 - 误 差 " 模 型 。 

TSE 模型 优点 在 于 可 以 分 离 出 状态 和 特质 两 
个 因子 ; ІНІМ, 它 将 往复 式 影响 精 确定 位 到 个 体 
内 层面 , 被 广泛 地 应 用 在 发 展 心理 学 和 临床 心理 
学 领域 中 。 如 抑郁 ， 可 分 为 特质 与 状态 两 个 层面 ， 
前 者 是 个 体 稳 定 的 特质 ， 随 时 间 推 移 保 持 稳定 ; 
而 后 者 是 情境 变量 ， 随 时 间 推 移 而 变化 (Hazel & 
Hankin, 2014; Kiken et al., 2015; Masselink et al., 
2018)。 分 离 出 稳定 的 特质 和 变化 的 状态 之 后 ,， 研 
究 者 可 以 根据 其 占 比 大 小 ,进行 适当 的 干预 ， 特 
别 是 从 特质 因素 的 根源 出 发 ， 寻 找 更 有 效 的 诊断 
方法 。 例 如 Masselink 等 人 (2018) 通 过 3 个 不 同样 
本 、 不 同时 间 跨 度 的 研究 发 现 , RI-CLM 建 模 指向 
“个 体内 仅 有 自尊 负 向 预测 抑郁 ”的 结论 ， 自尊 和 
抑郁 的 负 相关 是 源 于 特质 层面 ; ІП СІМ 建 模 则 
指向 “个 体内 自尊 和 抑郁 之 间 存 在 负 向 的 往复 式 
影响 ”的 结论 。 该 研究 的 差异 化 结果 恰好 证 实 了 贝 
克 认 知 理论 (Beck’s cognitive theory) 当 中 的 易 感 
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模型 (vulnerability model): 低 自 苯 的 个 体 容易 抑郁 ; 
而 本 来 就 难以 检测 的 创伤 模型 (scar model) 在 
RI-CLM 模型 中 被 证 伪 。 说 明 使 用 CLM 检测 的 灵 
敏 度 不 如 RI-CLM, 后 者 能 更 精准 地 定位 出 差异 
源 于 何 处 。 

1.4” 自 回归 模型 与 潜 增 长 模型 结合 


成 外 源 变量 ， 即 设 定 xi 和 yi 与 增长 因子 之 间 为 相关 
并 自由 估计 (Curran & Bollen, 2001; Jongerling & 
Hamaker, 2011; Zyphur, Voelkle, et al., 2020). Jb, 
如 果 将 斜率 项 px; 和 wwy; 限 定 为 0， 则 得 到 RI-CLM 
模型 (公式 7)， 故 RI-CLM 是 ALT 的 一 个 特例 。 该 
模型 需要 有 5 次 测量 才能 满足 时 变 参 数 模 型 的 识 


将 交叉 沾 后 模型 与 潜 增 长 模型 结合 ， 可 以 同 
时 考察 往复 式 影 响 与 增长 趋势 ， 其 建 模 方式 主要 
有 两 种 思路 。 第 一 种 建 模 方法 是 以 交叉 滞后 模型 
为 基础 模型 ， 在 其 基础 之 上 增加 潜 增 长 因子 。 其 
中 较 早 提出 来 的 是 Curran 和 Bollen (2001) 的 自 回 
JRA K i A (Autoregressive Latent Trajectory, 
ALT, Curran & Bollen, 2001)。 其 建 模 方式 如 下 ( 见 
图 6)。 

х = +E- Dihi + Paxi Уе) + Fault > D) 


Vir = My +E — Daxi + ra + уух) +d,,(t>1) 
在 结合 的 时 候 , 注意 到 xy 和 yi 定义 的 区 别 。 
因为 在 交叉 滞后 模型 当中 ,xi 和 yi 是 不 受 模型 内 
其 他 变量 影响 的 , 属于 “外 源 变量 ”; 而 在 定义 增 
长 因子 的 时 候 , x 和 入 被 当成 “指标 ”。 因 为 ALT 
的 载体 为 交叉 滞后 模型 ， 所 以 需要 将 xy 和 六 考虑 


别 要 求 ， 而 如 果 限 定 为 非 时 变 参数 ,4 次 测量 即 
可 。Usami 等 人 (2019) 也 把 ALT 模型 称 之 为 “累积 
增长 模型 ”增长 因子 被 看 成 是 两 个 “累积 因子 
(accumulating factors)”. 

ЯЕ SK A RS KER A ЖЕН, 
在 其 基础 之 上 增添 交叉 滞后 的 影响 系数 。 其 代表 
为 结构 化 残 差 潜 增 长 模型 (Latent Curve Model 
with Structured Residuals, LCM-SR, Bainter & 
Howard, 2016; Curran et al., 2014). 它 与 ALT 模型 
的 出 发 点 不 同 : ALT 出 发 点 是 自 回归 ， 以 此 添加 
个 体 间 变异 (从 个 体内 影响 中 分 离 出 个 体 间 变异 ); 
LCM-SR 出 发 点 是 潜 增 长 ， 以 此 添加 交叉 滞后 影 
响 ( 从 个 体 间 变异 中 分 离 出 个 体内 影响 )。 所 以 ， 
LCM-SR 不 纠结 于 第 一 次 测量 是 否 为 外 源 变量 ， 
直接 将 潜 变 量 因子 作用 于 所 有 测量 指标 ， 并 且 在 


图 6 自 回 归 洪 增长 模型 示意 图 
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测量 部 分 添加 交叉 滞后 效应 (图 7) 0 о 


Xit = xi + (7— D7 t P 1) FI GE не р + 


(10) 

Yi = My %(4-1) ты + Же 1) 十 JrE E yt dy 
其 中 ，s*xi ay All ety eye AE BER Far 
AR REMY 22201, AE AT] ZZ E Езу А БА 135 X Y 
后 参数 。 在 LCM-SR 中 , 更 清晰 地 看 到 观测 变量 的 
变异 被 分 解 成 了 被 试 间 和 被 试 内 两 部 分 。 被 试 间 
变异 表示 为 潜 增 长 模型 部 分 ,被 截 距 项 (mxai、7iv) 
和 和 斜率 项 (my2xi、 Nayi) ЕЖ; 被 试 内 变异 是 误差 项 
eža epf e* yi toD， 将 自 回归 和 交叉 清 后 的 影响 直 
接 定 义 在 测量 误差 之 间 (Bainter & Howard, 2016; 
Mulder & Hamaker, 2020)。 由 于 误差 项 本 身 也 是 洪 
变量 , 故 LCM-SR 模型 的 组 内 差异 建立 在 潜 变 量 上 。 
实际 中 如 何 定义 误差 之 间 的 关系 ? 如 果 将 e* 
看 成 是 Fx 的 一 个 特例 ， 那 e* 即 潜 变 量 因 子 , 定义 
其 观测 指标 x 和 yy 的 “测量 误差 ?为 0( 观 测 指 标 不 
再 含有 其 他 的 误差 因子 ),， 就 巧妙 地 将 理论 定义 操 
作 化 了 (参见 Mulder & Hamaker, 2020; Zyphur, 
Allison, et al., 2020)。LCM-SR 要 求 至 少 4 次 测量 


才能 满足 时 变 参 数 的 识别 要 求 ， 如 果 限 定 为 非 时 
变 参 数 则 需要 至 少 3 次 测量 。 

对 比 两 种 建 模 思路 可 知 ， 除 了 是 否 假设 第 一 
次 测量 为 外 源 变量 之 外 ，ALT 和 LCM-SR 的 建 模 
差异 集中 在 交叉 滞后 的 影响 是 哪 一 类 变量 所 指向 
的 。ALT 当中 的 交叉 滞后 影响 直接 建立 在 观测 变 

量 上 ， 即 不 考虑 测量 误差 ; 而 LCM-SR 当中 的 交 
叉 滞后 影响 建立 在 测量 误差 上 ， 即 潜 变 量 上 。 如 
果 将 ALT 模型 进行 拓展 ， 考 虑 其 每 个 观测 变量 的 
测量 模型 , 但 是 限定 这 些 测量 模型 的 “测量 误差 ” 
为 0， 即 可 得 到 LCM-SR 模型 ( 除 第 一 次 测量 为 外 
HEE) MAZ, LCM-SR 可 以 粗略 看 成 是 一 个 
限定 测量 误差 为 0 的 带 有 因 ен ALT 模型 。 

除 上 述 两 类 建 模 思想 以 外 ,， 潜 变化 分 数 模型 
(Latent Change Score ме, on 或 潜 差 异 分 数 
模型 (Latent Difference Score Model, LDS), 也 是 
同一 时 期 衍生 出 来 的 同时 关注 交叉 滞后 与 潜 增 长 
的 模型 (McArdle & Hamagami, 2001)。LCS 在 模型 
定义 上 , 首先 纳入 了 测量 误差 (公式 5), 然后 构造 


图 7 结构 化 残 差 潜 增长 模型 示意 图 
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“因素 差异 ”( 公 式 11), 并 且 用 差异 分 数 作 为 指标 
构造 增长 特质 (mx; ЖІ ту» 公式 12), ШІ 

Fa = Ёол) + AP uit 

Fun = Foray + APy 


yit yi(t-1 


ж 
AF in = AN xi + Buk 


xi(t-1 


(11) 


) + Yaf aai + d u(t > 1) 
AF i, = By N yi + By Аа) У.Р) + а, 2 1) 


模型 中 的 增长 特质 fii 和 ту 被 称 为 “斜率 ”， 
但 容易 将 上 述 公 式 进行 转换 ， 将 公式 12 带 入 公式 
11, 得 到 
Fi = а. ++ BOF at +у еу +d,,(t>1) 
Pi =Q hyi ++ By es ty ae ға,(4>1) 

可 知 , LCS (公式 13) 中 的 “斜率 因子 ”一 一 ii 
和 iy 一 一 其 实 本 质 上 仍然 是 ALT 当中 的 一 个 “ 截 
距 因 子 ”(Usami et al., 2015; Usami et al., 2016)。 
因子 载荷 а, 和 ay 代表 了 每 次 测量 的 权重 。 在 传 
统 LCS 模型 中 ,cs 和 ay, 被 设 定 为 非 时 变 的 ， 即 
ox=ay=1; 车 将 此 限定 放宽 ， 即 得 到 时 变 的 截 距 
项 Ох хі 和 ау туі» 也 被 称 为 三 重 变化 分 数 模型 
(triple change score, TCS, McArdle & Nesselroade, 
2014)。 这 样 就 可 以 在 LCS 基础 上 对 协 方差 进行 限 
Е, 得 到 FCLM 模型 (公式 6). ЖИ, ZARA K 
ERZ, 可 将 ЕСІМ 看 成 是 LCS 的 一 个 特例 
(Usami et al., 2015, 2016)。 除 此 之 外 ,二 者 的 自 回 
归 系 数 的 定义 也 略 有 差异 。 

同时 关注 潜 增 长 和 交叉 滞后 的 模型 在 模型 的 
解释 与 应 用 上 相对 复杂 ， 常 被 应 用 在 发 展 与 临床 
人 研究 当中 (Berry & Willoughby, 2017; Cole et al., 
2005; Curran et al., 2014; Ding et al., 2020; Malone 
et al., 2004)。 这 类 模型 解决 的 核心 问题 是 : 控制 
了 稳定 的 个 体 间 变异 以 后 , 个 体内 是 否 仍然 存在 
发 展 惰性 或 变量 之 间 的 往复 式 影 响 。 由 于 在 CLM 
中 往复 式 影 响 杂 糙 了 稳定 的 成 分 所 以 传统 模型 
的 参数 估计 与 整合 后 的 模型 结果 或 大 相 径 庭 ; 而 
如 果 考 察 的 目标 变量 确实 存在 这 种 “特质 ”( 研 究 
中 常 提 到 的 抑郁 、 焦 虑 、 自 苯 等 心理 学 成 分 ), 不 
经 控制 的 模型 会 高 估 这 种 影响 。 例 如 ，Berry 和 
Willoughby (2017) 的 研究 中 通过 两 组 真实 数据 对 
这 类 问题 进行 建 模 ， 探 索 儿 童 的 攻击 行为 与 父母 
的 体罚 ( 打 屁 股 ，spanking) 之 间 的 关系 。 发 现在 纳 


Ц 


CLM， 两 种 行为 之 间 存 在 交叉 滞后 影响 。 研 究 者 
指出 , 往复 式 影响 或 不 应 发 生 在 每 个 层面 。 因 为 
个 体 间 的 变异 应 是 稳定 的 非 时 变 的 ， 而 个 体内 变 
异 应 是 变化 的 、 波 动 的 ， 所 以 剔除 个 体 间 变异 对 
于 模型 的 可 靠 性 更 是 必要 的 前 提 (Berry & 
Willoughby, 2017)。 所 以 , 在 个 体 间 层面 ， 儿 童 的 
攻击 行为 受到 父母 体罚 行为 的 正 向 影响 ,然而 这 
并 不 会 发 生 在 个 体内 。 


2 ”模型 整合 与 拓展 


这 类 研究 的 核心 问题 在 于 如 何 分 解 变 异 。 洪 
增长 模型 所 提取 出 来 的 总 体 增长 趋势 被 看 成 是 个 
体 间 变异 ; 不 能 被 它 解释 的 ， 即 交叉 滞后 模型 探 
索 的 ， 是 个 体内 变异 (Bainter & Howard, 2016; 
Hamaker et al., 2015; Murayama et al., 2017). 在 提 
取 公 共 因 子 之 前 , 往复 式 影 响 发 生 在 个 体 的 所 
有 层面 ,但 是 当 分 离 出 特质 变量 之 后 ,往复 式 影 
响 被 定位 到 个 体内 层面 。 本 节 的 模型 整合 与 拓展 
思路 将 按照 追踪 研究 主要 解决 的 问题 来 展开 ， 
旨 在 解决 如 何 精确 定位 往复 式 影 响 和 增长 趋势 
的 问题 。 

首先 ， 研 究 者 若 关注 往复 式 影 响 , 可 以 建立 
交叉 滞后 模型 (CLM); 者 关注 个 体 增 长 趋势 ， 则 
可 建立 潜 增 长 模型 (LGM)。 然后, 分 别 以 CLM 或 
LGM 为 基础 模型 ， 添 加 相应 的 变异 分 解 成 分 。 诸 
如 测量 误差 ， 以 潜 变 量 指标 代替 观测 变量 进行 建 
模 , 或 者 根据 研究 需要 只 关注 截 距 因子 或 斜率 因 
子 。 大 考察 所 有 因素 ， 可 以 将 上 述 模型 进行 拓展 ， 
得 到 整合 的 模型 一 一 以 CLM 为 基础 模型 的 潜 变 量 
自 回归 潜 增 长 模型 (Latent Variable Autoregressive 
Latent Trajectory, LV-ALT, Bianconcini & Bollen, 
2018)， 以 及 以 LGM 为 基础 模型 的 因子 结构 化 潜 增 
长 模型 (Factor Latent Curve Model with Structured 
Reciprocals, FLCM-SR)。 

以 LV-ALT 为 例 (CLM 为 基础 模型 )， 其 模型 
建构 表示 如 公式 14 所 示 ( 见 图 8)。 

Fi =. + 0-Day + В.Р +7, Ее) жа,(>1) 


ж 


Қазу +(t—Drpy + ЖӘ Уи алу а, (4>) 
建 模 当中 的 各 个 参数 同 前 文 。 首 先 考 虑 测量 


с 


入 了 增长 模型 之 后 ， 体 罚 行 为 随 着 时 间 推 移 ， 频 
次 降低 (斜率 为 负 ), 但 是 攻击 行为 有 所 增加 ; 两 
种 行为 的 交叉 滞后 影响 不 显著 。 但 是 对 比 传统 


误差 (公式 5)， 自 回归 和 增长 曲线 均 在 真 分 数 Ға 
和 Fi 上 建 模 。 其 次 ,考虑 个 体 的 增长 形态 ， 即 截 
FE (1 All ту) AUR AE (oxi 和 туі), 其 系数 固定 为 特 


| 
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Ав 潜 变 量 自 回归 洪 增 长 模型 示意 图 


定 的 增长 形态 (其 与 第 一 次 测量 真 分 数 Ға 和 Fy 


量 误差 (公式 9)。 具 体操 作 层面 ,在 LCM-SR 基础 


设 定 为 相关 而 非 增长 形态 )。 在 此 基础 上 ， 建 立 自 
回归 影响 (86x, By) 和 往复 式 影响 (yx yy)。LV-ALT 
模型 识别 要 求 可 分 两 阶段 进行 : 结构 部 分 的 识别 
同 ALT， 时 变 效应 需要 至 少 5 次 测量 , 4 次 测量 要 
求 限定 线性 增长 ,3 次 测量 要 求 进一步 限定 自 回归 
系数 相等 ; 测量 部 分 的 识别 与 一 般 验 证 性 因素 分 
析 模 型 相同 (Bianconcini & Bollen, 2018). 

另 一 种 建 模 思 路 ,以 LGM 模型 为 基础 模型 
包含 了 所 有 分 解 变异 成 分 的 模型 即 因 子 结构 化 洪 
增长 模型 (FLCM-SR)。 其 构建 方法 如 公式 15( 见 图 
9)。 

Fai = Mau +0-0)т, + ВАР Я Р + dy (15) 
Pa = hyi +t- Dry + ВА + VF, eer жа, 

这 里 ,观测 指标 的 变异 首先 被 看 成 是 包含 真 
分 数 和 测量 误差 (公式 5), Манн a M Fyi 
的 变异 分 解 为 个 体 间 变 异 (包含 堆 距 mx; 和 yi、 
斜率 т 和 тоу). 个 体内 变异 (包含 Е р ЕЗ, 
4 的 自 回归 影响 和 往复 式 影 响 ) 和 残 差 (da 和 
dyit)。 注 意 到 公式 15 中 的 Fr 相当 于 公式 10 中 的 s*， 
这 意味 着 FCLM-SR 模型 中 将 结构 化 残 差 因 子 看 
成 是 真正 的 潜 变 量 ,追加 每 个 观测 指标 真正 的 测 


之 上 放宽 观测 指标 测量 误差 为 0 的 假设 即 可 (网 络 
版 附录 )。 

对 比 LV-ALT, 可 知 FCLM-SR 将 公式 14 中 
1 的 限定 放宽 , 忽略 第 一 次 测量 的 外 源 变量 性 
质 。 从 图 9 Ж, IE pixi 和 my: 直接 作用 于 хі 
和 уі 并 固定 载荷 为 1， 而 不 再 添加 因子 之 间 的 相 
关 ; 其 他 的 建 模 与 LV-ALT 相同 。 换言之 ， 两 类 模 
型 在 建 模 过 程 中 的 核心 问题 仍然 为 “是 否 将 第 一 
次 测量 看 成 外 源 变 量 ” 其 他 的 参数 设 定 都 只 是 
在 整体 框架 当中 的 特殊 案例 。 这 样 的 假设 在 实证 
当中 具有 可 操作 性 ， 可 以 让 第 一 次 观测 指标 从 属 
于 个 体 间 变异 ， 即 从 研究 一 开始 就 存在 固定 的 特 
质 。 因 为 实证 研究 的 首次 测量 点 并 不 一 定 是 特质 
发 展 当 中 真正 的 “原点 ” 或 存在 于 整个 研究 外 部 
的 真正 的 外 源 变量 。 故 可 在 实证 过 程 中 直接 用 第 
一 次 的 观测 指标 建 模 ， 由 此 一 来 对 于 结构 模型 识 
别 的 要 求 也 可 以 减少 一 次 (时 变 模型 4 次 测量 、 非 
时 变 模型 3 次 测量 )。 

Usami 等 人 (2019) 在 对 上 述 模型 进行 整合 的 
时 候 ， 将 变异 分 解 为 : (1) ll ee GR FE (measurement 
erronD)， 诸 如 因子 交叉 滞后 模型 、 特 质 -状态 -误差 模 
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型 中 对 测量 信 度 的 考察 ; (2) 成 分 方程 (decomposition 


思想 。 故 重新 对 多 变量 追踪 研究 的 系列 模型 所 


equation)， 即 潜 增 长 部 分 ,包括 只 有 随机 截 距 与 
同时 包含 截 距 和 和 斜率 的 情况 ; (3) 动 态 方程 
(dynamic equation)， 指 将 第 一 次 测量 看 成 外 源 变 
量 之 后 的 建 模 。 但 是 该 整合 框架 问题 在 于 ， 成 分 
方程 和 动态 方程 并 不 能 够 同时 建 模 一 一 实质 上 ， 
动态 方程 表示 了 以 交叉 滞后 模型 为 基础 的 建 模 思 
想 ， 而 成 分 方程 表示 了 以 潜 增 长 模型 为 基础 的 建 


H 

包含 的 成 分 进行 分 解 ， ОПА ЗА ie Je REA ЖЕ 
增长 模型 为 基础 模型 ， 直接 对 研究 者 都 比较 熟悉 
的 几 个 成 分 进行 增 减 建 模 。 拓 展 框架 如 表 1 所 示 。 

交叉 滞后 模型 当中 ,不 考虑 测量 误差 ， 也 不 
考虑 增长 形态 ， 建 模 最 简单 。 在 此 基础 上 ， 可 以 分 
别 考 察 测量 误差 (因子 交叉 滞后 模型 ) 和 随机 截 距 
(随机 截 距 交 叉 滞 后 模型 )， 或 将 二 者 同时 考虑 ( 特 


图 9 因子 结构 化 潜 增 长 模型 示意 图 


Жі 交叉 滞后 与 潜 增长 系列 模型 的 成 分 构成 一 览 


基础 模型 均值 


基础 模型 1: 交叉 滞后 模型 (CLMD 
因子 交叉 滞后 模型 (FE-CLMD) 

随机 截 距 交 叉 滞 后 模型 (RI-CLMD) 
特质 -状态 -误差 模型 (TSE/STARTS) 
自 回归 潜 增 长 模型 (ALT) 
潜 变 量 自 回 归 潜 增长 模型 (LV-ALT) = 


入 SoS 


基础 模型 2: 潜 增 长 模型 (LGMD) Е 
结构 化 残 差 潜 增 长 模型 (LCM-SR) Е 
因子 结构 化 潜 增 长 模型 (FLCM-SR) = 


WEAF 斜率 因子 滞后 残 差 测量 误差 
= у = 
- - м v 
v - v = 
v' - v v 
И А v - 
22 м v v 
v v = - 
v v v м 
v м v v 


注 .” 截 距 因 子 的 均值 为 0。 ”假设 第 一 次 观测 变量 为 外 源 变量 (与 增长 因子 为 相关 )。“ 限定 为 0。 
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质 -状态 -误差 模型 ); 在 此 基础 上 ， 加 入 增长 模型 学 期 和 八 年 级 春季 学 期 。 测 量 的 指标 有 直接 认 知 
( 自 回归 洪 增 长 模型 ) 及 其 测量 误差 ( 潜 变 量 自 回归 能 力 (包括 阅读 成 绩 、 数 学 成 绩 和 科学 成 绩 )、 间 接 
潜 增 长 模型 )。 男 一 个 思路 是 潜 增 长 模型 , 不 考虑 认 知 能 力 (包括 教师 对 学 生 认 知 能 力 的 整体 评价 ， 
交叉 滞后 ， 也 不 考虑 测量 误差 ， 建 模 也 相对 简单 。 采用 行为 量 表 计 分 )、 学 生 的 自我 描述 以 及 学 生 的 
在 此 基础 上 ,考虑 交叉 滞后 成 分 (结构 化 残 差 潜 增 自我 概念 和 控制 点 。 由 于 第 三 次 (一 年 级 下 期 ) 只 


长 模型 ), 或 进一步 考虑 测量 误差 (因子 结构 化 潜 增 
长 模型 )。 两 种 思路 最 终 统一 于 LV-ALT/FCLM-SR。 
研究 者 在 选用 模型 的 时 候 ， 应 该 以 某 一 类 模型 


对 其 中 30% 的 被 试 进行 了 调查 ,本 研究 中 没有 采 
用 这 一 次 的 调查 数据 ， 只 使 用 了 其 余 6 次 测量 的 
数据 ， 有 效 样本 21049 名 。 第 一 次 测试 时 学 生 的 


为 基础 模型 ， 在 此 基础 上 添加 相应 的 变异 分 解 平均 年 龄 68.47 Л, 标准 差 4.21 A, 最 大 值 79 
成 分 。 个 月 , 最 小 值 54 个 月 。 测 量 指标 采用 等 值 后 的 
ЕЕЕ 阅读 能 力 和 数学 能 力 IRT 分 数 ， 可 以 进行 跨 年 级 
3 RERA: 往复 式 影响 与 个 体 发 展 的 。 比较 。 测验 设计 进一步 的 信息 可 以 参考 相应 的 技 
模型 比较 、 
术 报 告 。 
3.1 ”数据 与 变量 数据 中 的 阅读 (1) 与 数学 (m) 能 力 的 不 同年 级 
采用 美国 国家 教育 统计 中 心 (National Center 的 均值 和 标准 差 以 及 相关 和 矩阵 如 表 2 所 示 。 
for Education Statistics, NCES) 提 供 的 早期 儿童 的 3.2 ”增长 шн 
追踪 研究 -幼儿 园 版 (Early Childhood Longitudinal 首先 ,探讨 增长 形态 。 根 据 前 人 研究 ( 刘 源 ， 
Study-Kindergarten cohort, ECLS-K, https://nces.ed. 刘 红 云 ， 2018)， 该 数据 的 增长 形态 或 存在 关键 转 
gov/ecls/)。 该 追踪 研究 调查 了 1998-1999 年 人 学 折 点 。 故 建立 三 个 模型 进行 比较 : 线性 增长 模型 、 
的 21409 名 学 生 , 用 长 达 9 年 的 追踪 调查 了 从 幼 二 次 增长 模型 和 多 阶段 增长 模型 。 将 能 力 增长 曲 
儿 园 到 八 年 级 的 学 生 的 智力 发 展 情况 和 来 自 不 同 线 和 原始 数据 进行 对 比 ， 如 图 10 所 示 。 从 结果 比 
层面 的 影响 因素 。 数 据 测 量 了 7 个 时 间 点 : 幼儿 较 中 可 以 看 出 ， 多 阶段 模型 对 于 数据 的 拟 合 效 果 
园 秋 季 学 期 ,幼儿园 春季 学 期 ， 一 年 级 秋季 学 期 ， 最 好 (网 络 版 附录 )， 故 考虑 多 阶段 模型 作为 潜 增 
一 年 级 春季 学 期 ， 三 年 级 春季 学 ， 五 年 级 春季 长 形态 。 
表 2 ECLS-K 学 生 阅 读 与 数学 能 力 发 展 均 值 、 标 准 差 与 相关 
变量 rl r2 r4 r5 r6 r7 ml m2 m4 m5 m6 m7 
rl 1.000 
r2 0.805 1.000 
r4 0.692 0.791 1.000 
r5 0.627 0.683 0.775 1.000 
r6 0.604 0.654 0.732 0.862 1.000 
17 0.551 0.569 0.625 0.761 0.802 1.000 
ml 0.784 0.735 0.689 0.673 0.657 0.597 1.000 
m2 0.703 0.765 0.708 0.690 0.675 0.608 0.834 1.000 
m4 0.598 0.658 0.729 0.691 0.677 0.606 0.7320 0.796 1.000 
m5 0.583 0.624 0.668 0.745 0.730 0.676 0.713 0.764 0.803 1.000 
m6 0.556 0.594 0.637 0.708 0.740 0.692 0.679 0.728 0.771 0.881 1.000 
m7 0.540 0.569 0.608 0.680 0.708 0.730 0.639 0.680 0.708 0.815 0.853 1.000 
M -1333 -0.754 0.093 0772 1.027 1.287 -1.181 -0.696 0.041 0.701 1.091 1.419 
SD 0.526 0.512 0.463 0.318 0.302 0.392 0.483 0.465 0.423 0.392 0.415 0.452 


TE. 相关 系数 均 达 到 p < 0.001 KF. 
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的 设 定 (网 络 版 附录 )。 


7 个 模型 的 模型 拟 合 与 参数 估计 结果 如 表 3 

ы 48 2... 
1 „#8 三 个 模型 可 以 看 成 是 逐步 排除 变异 的 过 程 ; 汪 
a б ОЧРАР ЫГ 41. 因子 并 排除 测量 误差 ， 模 型 拟 合 
Bs Б eo 升 ,然后 对 比 ALT 和 LV-ALT, enna 
за 一 阅读 二 次 х 阅读 实际 上 排除 了 测量 误差 , 模型 拟 合 有 所 提升 。 LCM-SR 
I -- 数学 二 次 о 数学 实际 和 FLCM-SR 的 对 比 和 上 一 对 类 似 , 后 者 在 排除 
Boe са a е BOG 了 测量 误差 之 后 ,模型 拟 合 得 到 大 幅度 提升 。 再 
年 级 次 , 对 比 添加 了 增长 趋势 的 4 个 模型 (FLCM-SR、 

图 10 ECLS-K 学 生 阅 读 、 数 学 能 力 增长 趋势 LCM-SR、LV-ALT ЖІ ALT), 整体 拟 合 都 不 如 没有 
添加 增长 趋势 的 3 个 模型 (TSE RI-CLM 和 CLM), 

зз ”往复 式 影响 的 模型 比较 这 也 和 增长 趋势 是 否 拟 合 数据 有 关 ( 网 络 版 附录 )。 
在 确定 了 增长 形态 的 基础 上 ， 添 加 交叉 滞后 最 后 ， 对 比 两 类 基础 模型 的 建 模 (是 否 包含 外 源 变 
部 分 。 建 立 如 表 2 所 示 的 7 个 模型 。RI-CLM 包 оң) 发 现 将 第 一 次 测量 考虑 成 内 源 变 量 的 模型 


含 一 个 随机 截 距 ,可 参照 LGM 建 模 的 方式 , 但 是 (FLCM-SR 和 LCM-SR) 要 比 考虑 成 外 源 变 量 的 模 
注意 RI-CLM 限定 潜 变 量 均值 为 0， 观 测 变量 截 距 型 (LV-ALT 和 ALT) 拟 合 更 优 。 从 实证 数据 中 综合 
自由 估计 (在 LGM 模型 中 一 般 限定 观测 变量 截 距 KE, ТЕМЕН Е ЖЕ, ЖЕЛІ 
为 0， 潜 变量 均值 自由 估计 )。 在 TSE 模型 中 ， 同 定 比较 复杂 , 参数 限定 较 多 ,整体 拟 合 也 较 差 。 整 理 
样 不 估计 截 距 因子 的 均值 ， 放 宽 指 标 截 距 自 由 估 来 看 , 特质 -状态 -误差 模型 最 适合 本 数据 。 

计 ; 同时 ,定义 测量 误差 .每 一 个 指标 都 定义 一 对 比 参数 估计 ,不 同 模型 得 出 一 致 性 结论 : 
个 淤 变量， 国定 载荷 为 1; 在 此 基础 上 ， 自 回归 和 自 回归 系数 比 交 叉 滞 后 系数 估计 值 大 。 阅 读 能 
交叉 滞后 模型 就 可 以 在 潜 变 量 上 进行 定义 了 。 的 自 回归 系数 比 数学 能 力 的 自 回归 系数 略 小 , 说 
ALT 考虑 多 阶段 增长 形态 (转折 点 在 第 三 次 测量 )、 明 数 学 能 力 更 倾向 于 惰性 发 展 ， 受 到 先前 状态 影 
自 回归 与 交叉 滞后 影响 : 注意 载荷 矩阵 A 只 影响 响 更 大 。 横 向 对 比 不 同 的 模型 ， 模 型 之 间 的 结果 
第 二 次 到 第 七 次 测量 ， 第 一 次 测量 与 潜 变 量 因子 差异 较 大 。 交 叉 清 后 参数 部 分 , RI-CLM 模型 得 出 
设 定 成 相关 并 自由 估计 。LV-ALT 可 将 TSE 的 语 了 负 向 影响 , 需 谨慎 下 结论 。 在 考虑 了 截 距 因 子 
名 与 ALT 的 语句 结合 即 可 (网 络 版 附录 )。 江 增长 滞后 ， 交 叉 滞 后 的 参数 估计 值 减 小 ; ПІЛІ Т 
模型 方面 ,LCM-SR 中 ,定义 测量 误差 为 潜 谈 量 并 率 之 后 ， 自 回归 部 分 略 有 增 大 。 

设 定 “误差 之 误差 ”为 0， 建 立 潜 变量 自 回 归 和 交 结合 本 研究 的 实际 情况 ,在 阅读 和 数学 的 增 
又 滞后 参数 。FLCM-SR 中 放宽 “误差 之 误差 "为 0 长 过 程 当中 ,存在 着 稳定 的 “能 力 特质 *， 即 从 幼 


表 3 交叉 滞后 潜 增 长 系列 模型 拟 合 (ECLS-K) 


拟 合 指数 FLCM-SR LCM-SR LV-ALT ALT TSE RI-CLM CLM 
AIC 43292.99 53602.70 61850.70 60321.31 17715.09 19356.16 19290.99 
BIC 43770.27 53976.57 62304.12 60806.54 18112.82 19777.75 19688.72 
卡 方 26373.52 37616.18 46737.20 45200.76 2258.98 4038.14 3834.60 
df 30 43 33 29 40 37 40 
CFI 0.841 0.774 0.719 0.704 0.987 0.976 0.975 
RSMEA 0.204 0.204 0.259 0.272 0.051 0.072 0.067 
SRMR 0.203 0.252 1.321 1.309 0.022 0.084 0.059 
注 : FLCM-SR: 因子 结构 化 潜 增 长 模型 ; LCM-SR: 结构 化 潜 增 长 模型 ; LV-ALT: 潜 变 量 自 回 归 潜 增长 模型 ; ALT: АЕА 


增长 模型 ; TSE: 特质 -状态 -误差 模型 ; RI-CLM: 随机 截 距 交 又 滞后 模型 ; СІМ: 交叉 滞后 模型 
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表 4 交叉 滞后 潜 增 长 系列 模型 参数 估计 (ECLS-K) 

we FLCM-SR LCM-SR LV-ALT ALT TSE RI-CLM CLM 
и Est. S.E. Est. S.E. Est. S.E. Ев. S.E. Est. S.E. Est. S.E. Est. S.E. 

自 回 归 交 叉 滞 后 影响 

R2 ON 

КІ 1.303 0.017 1.702 0.042 0.258 0.016 1.656 0.048 0.548 0.013 0.492 0.008 0.579 0.007 

МІ —0.391 0.013 -0.919 0.038 0.575 0.015 -0.655 0.041 0.349 0.022 0.031 0.010 0.286 0.008 

R4 ON 

R2 -1.423 0.042 0.531 0.018 -0.660 0.207 0.797 0.013 0.531 0.010 0.403 0.008 0.544 0.007 

M2 1.408 0.072 0.462 0.018 1.553 0.213 0.182 0.014 0.296 0.016 -0.003™ 0.009 0.250 0.007 

R5 ОМ 

R4 0.347 0.014 -0.670 0.065 -0.978 0.089 0.589 0.011 0.379 0.015 0.169 0.008 0.401 0.006 

M4 0.352 0.010 1.157 0.075 1.467 0.088 0.103 0.010 0.256 0.012 -0.021 0.007 0.203 0.006 

R6 ON 

RS 0.889 0.010 1.575 0.02 0.318 0.022 0.680 0.016 0.836 0.024 0.252 0.013 0.675 0.007 

MS 0.030 0.009 -0.466 0.014 0.526 0.020 0.071 0.011 0.045" 0.016 -0.061 0.009 0.150 0.006 

R7 ON 

R6 0.006" 0.029 1.745 0.004 -0.533 0.094 0.889 0.029 1.024 0.027 0.338 0.017 0.834 0.012 

M6 0.658 0.036 -0.400 0.015 1.493 0.082 0.171 0.019 0.107 0.022 0.034 0.011 0.206 0.009 

M2 ON 

МІ 0.274 0.013 1.401 0.041 0.633 0.018 0.186 0.058 0.920 0.017 0.290 0.010 0.708 0.007 

КІ 1.200 0.049 -0.592 0.043 0.186 0.018 0.644 0.063 -0.042 0.015 0.102 0.007 0.113 0.006 

M4 ON 

M2 1.046 0.041 0.497 0.017 1.977 0.185 0.573 0.018 0.832 0.025 0.203 0.009 0.647 0.007 

R2 -0.952 0.027 0.342 0.013 -1.063 0.178 0.322 0.018 -0.022" 0.019 0.043 0.007 0.096 0.006 

MS ON 

M4 0.409 0.011 -0.396 0.037 -0.557 0.095 0.777 0.017 0.806 0.033 0.197 0.008 0.628 0.008 

R4 0.374 0.020 0.661 0.031 1.252 0.096 0.114 0.015 0.045" 0.022 0.010" 0.007 0.150 0.007 

M6 ON 

MS -0.049% 0.017 0.905 0.009 0.961 0.020 1.043 0.016 1.105 0.020 0.373 0.010 0.840 0.007 

RS 1.221 0.023 0.101 0.007 -0.008 0.024 -0.161 0.016 -0.112 0.021 -0.106 0.012 0.149 0.009 

M7 ON 

M6 0.726 0.017 0.884 0.023 1.735 0.074 1.053 0.024 0.970 0.013 0.363 0.011 0.792 0.009 

R6 0.0297 0.009 —0.005" 0.003 -0.507 0.089 -0.077 0.026 0.073 0.019 -0.046 0.015 0.254 0.012 

限制 参数 影响 * 

Bx 0.459 0.019 0.161 0.026 0.371 0.024 0.578 0.007 0.419 0.011 0.322 0.007 0.536 0.003 

Үх 0.579 0.023 0.349 0.014 0.339 0.024 0.074 0.009 0.340 0.018 -0.072 0.006 0.235 0.003 

By 0.210 0.016 0.763 0.034 0.501 0.015 0.564 0.013 0.930 0.014 0.214 0.006 0.721 0.003 

Yy 0.434 0.016 0.033 0.037 0.196 0.016 0.177 0.009 -0.017 0.011 -0.044 0.005 0.116 0.003 

注 . ”限制 参数 影响 表示 将 对 应 自 回归 系数 、 交 叉 滞后 系数 限定 为 相等 ,估计 平均 效应 ， 即 非 时 变 模型 。 


除 标注 以外, 均 达 到 p < 0.001 KF. 
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儿 园 到 初中 的 发 展 中 保持 个 体 间 的 稳定 性 ,在 此 
基础 之 上 , 个 体内 发 展 惰性 主要 存在 于 低 年 级 ; 
在 高 年 级 时 个 体内 的 澡 留 影响 较 小 ， 容 易 受 到 其 
他 能 力 的 影响 ; 阅读 和 数学 能 力 之 间 也 具有 往复 
式 影响 。 从 增长 形态 来 说 ， 多 阶段 增长 模型 虽然 
比 线性 或 二 次 模型 有 较 好 的 拟 合 , 但 是 不 及 单一 
随机 截 距 因子 模型 拟 合 优良 。 综 合 看 来 ， 可 采用 
TSE 来 建 模 。 


4 讨论 与 总 结 


41 BEBE: 如 何 分 解 变异 

本 文 探 讨 了 多 元 追踪 研究 中 的 不 同 模型 ， 对 
模型 的 构建 和 应 用 进行 了 梳理 ， 并 且 将 模型 进行 
整合 ， 拓 展 出 具有 操作 化 意义 的 模型 框架 。 结 合 
以 往 研究 ，Usami 等 人 (2019) 的 框架 并 没有 考虑 
AILT 模 型 的 测量 部 分 ， 以 及 LCM-SR 固定 “误差 之 
误差 ?为 0 的 假设 .虽然 Bianconcini 和 Bolle (2018) 
的 模型 考虑 了 ALT 的 测量 误差 , 但 是 目前 的 研究 
大 多 以 CLM 为 基础 模型 进行 建 模 ,忽略 了 LGM 
为 基础 模型 的 建 模 和 二 者 之 间 的 转换 。 故 本 文 综 
合 梳理 了 两 个 建 模 思 路 ,各 自 增 减 相 应 的 变异 成 
分 , 得 到 不 同 的 模型 。 

以 本 文 整合 框架 中 的 FLCM-SR Wil, 首先 
排除 个 体 间 变 异 ， 即 潜 增 长 部 分 。 这 部 分 描述 了 
个 体 随时 间 发 展 的 稳定 特征 ， 既 描述 了 初始 特质 
GER), 又 刻画 了 增长 规律 (成 熟 )。 建 模 过 程 中 ， 
潜 增 长 部 分 可 以 被 看 成 是 一 个 需要 被 控制 的 特征 ， 
由 于 个 体 “ 基 线 ”" 和 “成 熟 ” 的 存在 ,排除 掉 一 部 分 
自然 发 展 的 规律 是 合理 的 (Berry & Willoughby, 


归 即 滞留 影响 或 个 体 惰性 , 研究 者 应 将 其 和 前 一 
段 提 到 的 个 体 间 的 “自然 成 熟 " 进 行 区 分 (Berry & 
Willoughby, 2017; Curran & Bauer, 2011; Hamaker 
et al., 2015; Voelkle et al., 2014). “滞留 影响 ”是 指 
个 体内 当前 的 状态 多 大 程度 上 取决 于 他 /她 的 已 
有 状态 ,， 即 变化 的 顽固 程度 ， 是 个 体 不 变 的 惯性 ; 
而 “自然 成 熟 ?是 个 体 间 的 发 展 规律 ， 就 像 诸 多 心 
理学 教材 当中 书写 的 “既定 内 容 ” 一 样 。 以 本 文 实 
证 研究 的 阅读 发 展 为 例 ,“ 自 然 成 熟 " 即 潜 增 长 部 
分 是 一 般 的 阅读 发 展 规律 一 一 幼儿 园 到 一 年 级 ， 
阅读 发 展 迅速 ,进入 小 学 之 后 ， 阅 读 发 展 速度 降 
下 来 ,而 “ 灌 留 影响 ”是 指 ,诸如 一 年 级 的 阅读 能 
受到 幼儿 园 阅 读 能 力 的 影响 ， 八 年 级 的 阅读 能 力 
受到 五 年 级 阅读 能 力 的 影响 。 这 种 滞留 影响 在 一 
些 灵活 的 SEM 结构 中 可 以 自由 估计 ， 比 如 一 年 级 
阅读 能 力 受 幼儿 园 阅 读 能 力 的 影响 较 大 ， 而 八 年 
级 阅读 能 力 受 五 年 级 阅读 能 力 影 响 较 小 。 所 以 ， 
此 类 模型 区 分 了 个 体内 和 个 体 间 效应 , 解决 了 发 
展 当 中 普遍 规律 与 阶段 差异 的 问题 , 更 强调 个 体 
在 环境 当中 的 “可 塑性 ”。 

下 面 一 部 分 是 交叉 滞后 影响 ， 即 往复 式 影 响 ， 
这 部 分 结果 或 许 是 研究 者 最 为 关注 的 焦点 之 一 。 
在 排除 了 个 体 间 变异 之 后 ， 交 叉 滞后 的 影响 或 被 
淹没。 原因 有 可 能 是 在 比较 传统 的 模型 中 ， 交 叉 
滞后 的 影响 杂 熔 了 发 展 、 个 体 间 增 长 规律 等 因素 ， 
使 得 交叉 滞后 的 影响 被 高 佑 (Berry & Willoughby, 
2017)。 但 是 , 这 种 “高 估 ” 并 非 完 全 不 可 接受 ， 
为 这 种 杂 灶 了 其 他 因素 的 影响 是 在 某 个 水 平 层 面 
上 的 解读 ， 可 以 把 它 看 成 是 “整体 影响 "此 类 研究 


2017)。 特 别 地 , 在 建立 潜 增 长 模型 的 时 候 ， 如 果 
测量 的 变量 之 间 没 有 经 过 等 值 , 则 斜率 因子 的 意 
义 不 明 确 , 或 只 能 反映 个 体 在 群体 当中 的 相对 变 
化 (Liu et al., 2020)。 在 本 文 实 证 应 用 的 建 模 过 程 
т, 除了 CLM 之 外 ,其实 都 考虑 了 个 体 间 的 成 熟 
ALA, AAW RI-CLM 和 TSE RIIA TRIE, 
ALT/LCM-SR 和 LV-ALT/FLCM-SR 同时 加 入 了 截 
距 和 斜率。 本 例 中 加 入 斜率 的 模型 估计 结果 稍 过 
于 只 加 截 距 的 模型 ， 可 能 是 多 阶段 形态 对 本 组 数 
据 的 拟 合 并 非 最 优 的 结果 而 导致 。 

在 排除 了 成 熟 影响 之 后 ,就 可 以 考察 个 体内 
的 影响 了 ， 即 FLCM-SR 中 的 “SR (结构 化 往复 式 
影响 )”。 在 这 里 ,个 体内 的 影响 主要 由 自 回归 和 交 
又 滞后 参数 两 部 分 构成 。 特 别 地 ,个 体内 的 自 回 


= 


倾向 于 将 往复 式 影响 定位 到 个 体内 层面 , 或 可 以 
更 精确 地 表达 如 本 文 一 开始 就 定义 的 “变化 是 如 
何 引起 变化 ”的 含义 。 

最 后 的 议题 是 是 否 考 虑 测量 误差 ， 即 每 次 观 
测 变量 的 测验 信和 度 。 本 例 中 LV-ALT 和 FLCM-SR 
的 拟 合 分 别 略 优 于 ALT 和 LCM-SR, 可 见 考 虑 了 
测量 信 度 会 一 定 程度 上 提升 模型 估计 质量 。 值 得 
注意 的 是 ,本 例 中 每 次 测量 都 是 单一 指标 ， 所 以 
测量 信和 度 是 一 个 固定 值 。 如 果 每 次 的 测量 有 多 个 
指标 ， 或 者 达到 模型 识别 的 要 求 , 测量 误差 才 真 
正 被 模型 估计 。 可 以 参考 Bianconcini 和 Bollen 
(2018) 提 出 的 LV-ALT 或 Mulder 和 Hamaker (2020) 
提出 的 多 指标 RI-CLM 框架 ,建构 多 指标 测量 模 
型 。 综 合 看 来 ， 如 果 有 人 先 验 信息 ， 加 入 测量 信和 度 的 


еш 
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模型 拟 合 会 更 精确 ,但 是 在 只 有 单一 指标 的 情况 
下 ， 转 于 模型 识别 等 因素 ， 必 须 人 为 限定 测量 信 
度 。 此 外 , 在 多 指标 跨 时 测量 时 ， 跨 时 测量 等 价 性 
也 是 需要 研究 者 提前 考察 的 议题 (Masselink et al., 
2018)。 
4.2 ”模型 选择 的 议题 

上 述 论证 的 变异 分 解 理论 建构 应 该 是 研究 者 
考虑 的 模型 选择 的 重要 依据 之 一 。 尽 管 如 此 ， 模 
型 选择 的 过 程 中 , 确 有 一 些 不 稳定 因素 会 影响 研 
究 者 的 使 用 。 比如, 研究 者 并 不 知道 “ 真 ”模型 ， 即 
理论 上 真正 的 变异 应 该 如 何 分 解 ， 所 以 无 法 选择 
“正确 ”的 模型 ,事实 上 , 一 些 模 拟 研究 的 结果 也 显 
ж, 通过 不 同 定 义 框架 来 生成 的 数据 ， 采用 “ 正 
确 ” 模 型 的 估计 反 真 率 也 几乎 只 有 六 七 成 (Usami 
et al., 2015; 2016)。 这 也 表明 ,在 对 上 述 模型 的 选 
择 上 ,“ 正 确 ” 模 型 能 被 选 出 来 的 概率 并 不 是 太 高 ; 
而 “错误 ”模型 使 用 之 后 ， 其 实 对 结果 解释 并 不 会 
造成 太 大 的 影响 。 所以， 在 实证 研究 中 ， 并 不 是 模 
型 越 复杂 拟 合 越 佳 ， 系 列 的 拟 合 指标 是 研究 者 选 
择 模型 的 首要 参考 依据 (Usami et al., 2015). 

建立 复杂 的 模型 容易 , 但 是 选择 合适 的 模型 

不 容易 。 人 研究 者 应 首先 根据 理论 基础 设 定 合理 
的 模型 。 如 果 某 个 特质 因素 的 确 是 稳定 特质 ， 则 
考虑 包含 随机 截 距 的 模型 以 提出 稳定 特质 ; 如 果 
测量 时 间 较 长 且 确 实 有 发 展 ( 如 能 力 )， 再 考虑 增 
长 因子 。 一 般 而 言 ， 提 取出 稳定 特质 的 模型 会 具 
有 较 好 的 模型 拟 合 ， 而 增长 因子 与 交叉 滞后 部 分 
夹杂 在 一 起 会 带 给 研究 者 一 些 解 释 困 难 。 其 次 ， 
通过 备 选 模型 之 间 的 拟 合 指标 比较 来 选择 模型 是 
应 用 研究 者 应 该 参考 的 主要 评判 标准 之 一 。 此 类 
模型 由 于 复杂 程度 差异 较 大 ， 前 人 模拟 研究 结果 
建议 更 多 采用 CFI 和 RMSEA 来 选择 模型 ; 信息 
函数 的 表现 并 不 理想 , Н BIC 容易 惩罚 复杂 模型 
而 选择 更 简单 的 模型 (Usami et al., 2015)。 综 合 以 
往 人 研究 的 建议 (Berry & Willoughby, 2017; Usami 
et al.，2019) 和 本 文 实证 数据 的 结果 , 应 用 者 应 使 
用 诸如 CFI, TLI, RMSEA, SRMR 以 及 AIC 等 
指标 进行 综合 判断 。 
42 ”其 他 议题 与 研究 展望 

对 于 模型 的 选择 ， 近 期 有 学 者 提出 “结构 方 
程 模 型 树 (SEM tree)*” 和 “结构 方程 模型 森林 (SEM 
forest)”。 借 鉴 了 管理 心理 学 当中 的 决策 树 概念 ， 
来 对 模型 的 选择 提供 证 据 (Brandmaier et al., 2016; 


Brandmaier et al., 2013)。 未 来 可 以 通过 模拟 研究 
来 给 出 模型 选择 的 建议 ， 比 如 如 何 选 择 真 模型 ， 
模型 误 设 的 代价 ， 以 及 个 体 间 /个 体内 差异 达到 何 
种 程度 的 时 候 推 荐 哪 类 模型 等 。 

另 一 方面 ,时 变 模型 的 构建 可 以 使 得 模型 估 
计 更 为 精确 。 比 如 自 回归 模型 中 的 非 时 变 影响 在 
SEM 框架 下 可 以 将 参数 进行 一 定 程度 的 放宽 ; И 
留影 响 B, 在 不 同 测量 次 数 之 间 可 以 自由 估计 , 模 
型 限定 更 少 ， 可 以 有 效 提高 模型 整体 拟 合 。 
Bringmann 等 人 (2018) 提 出 了 时 变 自 回归 模型 
(time-variant autoregressive model), 假设 自 回 归 
影响 在 不 同 次 的 时 间 间 隔 不 同 (Bringmann et al., 
2018)。 若 设 定 清 留 效应 的 非 时 变性 的 理论 假设 并 
不 确定 ， 则 可 以 放宽 该 假设 ,整体 拟 合 更 优 (Ding 
et al., 2020)。 但 是 需要 注意 参数 识别 的 问题 。 

此 外 ,研究 者 在 追踪 研究 中 根据 不 同 测量 时 
间 间 隔 ， 对 于 时 间 变 量 的 精度 要 求 会 有 所 不 同 。 
在 SEM 框架 下 ,载荷 矩阵 A 大 多 都 固定 测量 时 
间 。 这 会 增 大 估计 误差 ， 降 低 模型 的 敏感 度 (Liu 
et al., 2015; Sterba, 2014)。 在 此 基础 上 ,借鉴 多 水 
平 模型 框架 ， 可 以 放宽 测量 时 间 固 定 的 假设 ， 
SEM 也 逐渐 开始 使 用 连续 时 间 模 型 (continuous 
time models，CTMD)， 将 被 试 参 与 的 时 间 定 义 为 一 
个 随机 数 ， 而 非 固 定 的 测量 次 数 (Deboeck & 
Preacher, 2016; Hamaker & Wichers, 2017; 
Schuurman et al., 2016)。 在 新 版 本 的 SEM 软件 (如 
Mplus 8.0) 中 ， 这 类 问题 也 已 经 得 到 技术 支撑 ， 可 
以 供 研究 者 灵活 选用 。 
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A unification and extension on the multivariate longitudinal models: 
Examining reciprocal effect and growth trajectory 


LIU Yuan 
(School of Psychology, Southwest University; Key Laboratory of Cognition and Personality 


(Southwest University), Ministry of Education, Chongging 400715, China) 


Abstract: When conducting the multivariate longitudinal studies, reciprocal relationship and latent 
trajectory are two of the focusing issues. These two issues could be flexibly combined by other research 
questions, such as the measurement error, the random factor, as well as the combination of the above topics. 
Such a combination yields a more complex model definition exploring the longitudinal relations, such as 
factor cross-lagged model, random-intercept cross-lagged model, trait-state-error model, autoregressive 
trajectory model, etc. In the study, a factor latent curve model with structured reciprocals model was built as 
an extension and unified framework including all the components discussed above. The empirical dataset, 
Early Childhood Longitudinal Survey-Kindergarten (ECLS-K), was used as an illustrating example. Results 
indicated that the trait-state-error model best described the data. Finally, we summarized how the results 
could be interpreted and offered suggestions on model selection for the researchers. 
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附录 1: 线性 、 二 次 、 多 阶段 增长 模型 结果 参数 估计 


线性 增长 模型 (LGM) 二 次 增长 模型 (QGM) 多 阶段 增长 模型 (PGM) 

Est. S.E. Est. S.E. Est. S.E. 
因子 均值 
INTR —0.889 *** 0.004 ПТТ жж 0.004 -1.321%%% 0.004 
SLPR 0.353 *** 0.001 0.729 *** 0.002 1.158 *** 0.003 
QUAR —0.055 *** <0.001 
SLPR2 —1.009 *** 0.004 
INTM —0.865 *** 0.003 —1.024 *** 0.003 —1.167 *** 0.003 
SLPM 0.363 *** 0.001 0.640 *** 0.001 0.914 *** 0.002 
QUAM —0.042 *** <0.001 
SLPM2 —0.703 *** 0.003 
因子 方差 
INTR 0.205 *** 0.003 0.244 *** 0.003 0.240 *** 0.003 
SLPR 0.004 *** <0.001 0.015 *** 0.001 0.026 *** 0.001 
QUAR <0.001 *** <0.001 
SLPR2 0.014 *** 0.002 
INTM 0.184 *** 0.002 0.196 *** 0.002 0.200 *** 0.002 
SLPM 0.003 *** <0.001 0.008 *** <0.001 0.008 *** 0.001 
QUAM <0.001 *** <0.001 
SLPM2 —0.003 *** 0.001 
测量 误差 
rl 0.288 *** 0.004 0.131 *** 0.002 0.045 *** 0.002 
12 0.065 *** 0.002 0.041 *** 0.001 0.059 *** 0.001 
т4 0.266 *** 0.003 0.120 *** 0.002 0.180 *** 0.003 
15 0.202 *** 0.003 0.014 *** 0.001 0.026 *** 0.001 
гб —0.009 *** 0.001 0.069 *** 0.001 0.008 *** <0.001 
17 0.657 *** 0.010 0.180 *** 0.008 0.090 *** 0.002 
ml 0.163 *** 0.002 0.083 *** 0.001 0.042 *** 0.001 
m2 0.035 *** 0.001 0.029 *** 0.001 0.035 *** 0.001 
m4 0.179 *** 0.002 0.088 *** 0.001 0.082 *** 0.001 
m5 0.120*** 0.002 0.017 *** 0.001 0.030 *** 0.001 
m6 0.002 0.001 0.041 *** 0.001 0.024 *** 0.001 
m7 0.617 *** 0.010 0.045 *** 0.004 0.134 *** 0.003 
因子 相关 
SLPR WITH 
INTR -0.020%%% <0.001 -0.052%%% 0.001 —0.069 *** 0.002 
QUAR 0.004 *** <0.001 
SLPR2 0.060 *** 0.002 
SLPR WITH 
QUAR —0.001 *** <0.001 
SLPR2 —0.020 *** 0.001 
SLPM WITH 
SLPM —0.010*** <0.001 —0.023 *** 0.001 —0.031 *** 0.001 
QUAM 0.002 *** <0.001 
SLPM2 0.026 *** 0.001 
SLPM WITH 
QUAM —0.001 *** <0.001 
SLPM2 —0.003 *** 0.001 


iE: *** p< 0.001 
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附录 2: 潜 变 量 自 回 归 潜 增长 模型 Mpls 语句 与 示意 


TITLE: this is an example of LV-ALT 
DATA: FILE IS eclsk thetatalt.dat; 
VARIABLE: NAMES ARE 
ClR4RTHT К СІКАМТНТ К CIRGTHT_R C2R4RTHT_R C2R4MTHT_R C2RGTHT_R C3R4RTHT_ R 
C3R4MTHT_R C3RGTHT_R C4R4RTHT_R C4R4MTHT_R C4RGTHT_R С5КАКТНТ К CSR4MTHT К 
CSR2STHT К C6R4RTHT_R C6R4MTHT_R C6R2STHT_R С7КАКТНТ К С7КАМТНТ К C7R2STHT_R 
PILEARN P2LEARN P4LEARN TILEARN T2LEARN T4LEARN TS5LEARN T6LEARN; 
USEVARIABLES ARE 
C1R4RTHT_R C2R4RTHT_R C4R4RTHT_R С5КАКТНТ К C6R4RTHT_R C7R4RTHT_R 
СІКАМТНТ К C2R4MTHT_R C4R4MTHT_R С5КАМТНТ К СӨКАМТНТ К С7КАМТНТ К; 
MISSING АКЕ ALL (-9); 
MISSING ARE ALL (-8); 
MISSING ARE ALL (-7); 
MISSING ARE ALL (-1); 
ANALYSIS: MODEL = NOCOVARIANCES; 
estimator=mlr; 
miteration = 5000;convergence=0.01; 
MODEL: 
!define the growth (排除 第 一 次 测量 为 外 源 变 量 后 , 定义 多 阶段 潜 增 长 模型 ， 前 三 次 测量 为 斜率 1， 
和 三 次 以 后 为 斜率 2); 
INTR by C2R4RTHT_R@1 C4R4RTHT_R@1 С5КАКТНТ К(@1 COR4RTHT_R@1 C7R4RTHT R@!1; 
SLPR by C2R4RTHT_R@0.5 САКАКТНТ Е(а)1.5 
CS5R4RTHT_R@3.5 C6R4RTHT_R@S.5 СТКАКТНТ К(а)8.5; 
SLPR2 by C2R4RTHT_R@0 САКАКТНТ R@O 
CS5R4RTHT_R@2 C6R4RTHT_R@4 C7R4RTHT_R@7; 
INTM by C2R4MTHT_R@I1 САКАМТНТ R@1 С5КАМТНТ R@1 СӨКАМТНТ R@1 C7R4MTHT К(а)1; 
SLPM by С2КАМТНТ К(@0.5 САКАМТНТ К(@1.5 
CSR4MTHT_R@3.5 СӨКАМТНТ К(а)5.5 С7КАМТНТ К(а)8.5; 
SLPM2 by C2R4MTHT_R@0 САКАМТНТ К(а)0 
CSR4MTHT_R@2 C6R4MTHT_R@4 C7R4MTHT_R@7; 
[C2R4RTHT_R-C7R4MTHT_R@0]; 
INTR SLPR INTM SLPM slpr2 slpm2; 
[INTR SLPR INTM SLPM slpr2 slpm2]; 
intr with slpr slpr2; slpr with slpr2; 


ne 


intm with slpm slpm2; slpm with slpm2; 

!define state factor (定义 状态 变量 ， 即 每 次 测量 的 信和 度 ); 

R1 BY СІВАКТНТ К(@1; R2 BY C2R4RTHT R@1; R4 BY C4R4RTHT R@I; 

R5 BY C5R4RTHT_R@1; R6 BY C6R4RTHT_ R@1; R7 BY C7R4RTHT R@I; 
МІ BY CIR4MTHT_ R@1; M2 BY C2R4MTHT R@1; M4 BY C4R4MTHT R@!1; 
М5 BY CS5R4MTHT_R@1; M6 BY C6R4MTHT_ R@1; M7 BY СТКАМТНТ R@I1; 
!define the exogenous indicators (定义 第 一 次 测量 的 因子 为 外 源 变量 ); 
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КІ with intr slpr slpr2; 

МІ with intm slpm slpm2; 

!define measurement error (定义 误差 变量 ); 
C1R4RTHT_R-C7R4RTHT К (1); 
C1IR4MTHT_R-C7R4MTHT К (2); 

!define auto-regression and cross-lagged regression (ХЕ Ў А EIJI- 48 Shir a2 Md); 
R2 ON ВІ; R4 ON В2; R5 ОМ R4; 

R6 ОМ RS; R7 ON R6; 

M2 ON МІ; M40NM2; M5ON M4; 

М6 ОМ М5; M7ON M56: 

M2ONRI: М40МЕ2; MS5SONR4; 

М6 ОМ А5; М70МЕб; 

R2 ONM1; Р40ОММ2; В5ОММд4; 
В6ОММ5; R7ON M6; 

!define the residuals (定义 外 源 变量 方差 、 被 试 内 残 差 ); 
R1; R2-R7(3); 

М1; М2-М7(4); 

!define the correlation between residuals (定义 因子 之 间 协 方差 /相关 ); 
КІ WITH МІ; R2 WITH М2; R4 WITH M4; 

RS WITH М5; R6 WITH М6; R7 WITH M7; 

OUTPUT: STDYX;tech4;tech1; 
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附录 з: 因子 结构 化 潜 增长 模型 Mplus 语句 与 示意 图 


TITLE: this is an example of FLCM-SR 
DATA: FILE IS eclsk_theta+alt.dat; 
VARIABLE: NAMES ARE 
ClR4RTHT К СІКАМТНТ К CIRGTHT_R C2R4RTHT_R C2R4MTHT_R C2RGTHT_R C3R4RTHT R 
C3R4MTHT_R C3RGTHT_R C4R4RTHT_R C4R4MTHT_R C4RGTHT_R CSR4RTHT К С5КАМТНТ К 
CS5R2STHT_R C6R4RTHT_R C6R4MTHT_R C6R2STHT_R C7R4RTHT_R C7R4MTHT_R 
C7R2STHT_R 
PILEARN P2LEARN P4LEARN TILEARN T2LEARN T4LEARN TS5LEARN T6LEARN; 
USEVARIABLES ARE 
СІКАКТНТ К C2R4RTHT_R C4R4RTHT_R С5КАКТНТ К C6R4RTHT_R C7R4RTHT_R 
СІКАМТНТ К C2R4MTHT_R C4R4MTHT_R С5КАМТНТ К C6R4MTHT_R С7КАМТНТ К; 
MISSING АКЕ ALL (-9); 
MISSING ARE ALL (-8); 
MISSING ARE ALL (-7); 
MISSING ARE ALL (-1); 
ANALYSIS: MODEL = NOCOVARIANCES; 
estimator=mlr; 
miteration = 5000;convergence=0.01; 
MODEL: 
!define the growth (定义 多 阶段 法 增长 模型 ， 前 三 次 测量 为 斜率 1, 第 三 次 以 后 为 斜率 2); 
INTR by CIR4RTHT_R@1 C2R4RTHT_R@1 САКАКТНТ К@ 1 С5КАКТНТ К(а)1 C6R4RTHT R@1 
C7R4RTHT К@1; 
SLPR Бу СІКАКТНТ К(@0 C2R4RTHT_R@0.5 C4R4RTHT К(а)1.5 
CS5R4RTHT_R@3.5 C6R4RTHT_R@S.5 СТКАКТНТ К(а)8.5; 
SLPR2 by CIR4RTHT_R@0 C2R4RTHT_R@0 C4R4RTHT_R@O 
С5КАКТНТ К(а)2 C6R4RTHT_R@4 C7R4RTHT_R@7; 
INTM by СІКАМТНТ R@1 C2R4MTHT R@1 C4R4MTHT R@1 С5КАМТНТ R@1 C6R4MTHT R@1 
C7R4MTHT_R@1; 
SLPM by СІКАМТНТ К@0 C2R4MTHT_R@0.5 C4R4MTHT_R@1.5 
CSR4MTHT_R@3.5 СӨКАМТНТ К(а)5.5 СТЕАМТНТ К(@8.5; 
SLPM2 Бу CIR4MTHT_R@0 C2R4MTHT_R@0 САКАМТНТ К(а)0 
CSR4MTHT_R@2 C6R4MTHT R@4 C7R4MTHT_R@7; 
[CIR4RTHT_R-C7R4MTHT_R@0]; 
INTR SLPR INTM SLPM slpr2 slpm2; 
[INTR SLPR INTM SLPM slpr2 slpm2]; 
intr with slpr slpr2; slpr with slpr2; 


ішіп with slpm slpm2; slpm with slpm2; 

!define state factor (定义 状态 变量 ， 即 每 次 测量 的 信和 度 ); 

КІ BY CIR4RTHT_R@1; R2 BY C2R4RTHT_R@1; КА BY C4R4RTHT R@!1; 
R5 BY С5КАКТНТ К@ 1; R6 BY C6R4RTHT_R@I1; R7 BY C7R4RTHT R@!1; 
МІ BY СІКАМТНТ К@1; M2 BY C2R4MTHT_R@|1; M4 BY C4R4MTHT R@I1; 
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М5 BY CS5R4MTHT_R@1; M6 BY C6R4MTHT R@!1; M7 BY СТКАМТНТ R@!1; 
!define measurement error (定义 误差 变量 ); 

СІКАКТНТ К-С7КАКТНТ К; 

СІКАМТНТ К-С7КАМТНТ К; 

!define auto-regression and cross-lagged regression (ХЕ Ў А 1919-5725 іг ЛЕ ЖЕРІН); 
R2 ON R1; R4 ON R2; R5 ON R4; 

R6 ON R5; R7 ON R6; 

M2 ON M1; M4 ON M2; M5 ON M4; 

M6 ON M5; M7 ON M6; 

M2 ON R1; M4 ON R2; M5 ON R4; 

M6 ON R5; M7 ON R6; 

R2 ON МІ; R4 ON M2; RS ON M4; 

R6 ON M5; R7 ON M6; 

Idefine the residuals (72 X WIA A FEF); 

R1; R2-R7(3); 

МІ; М2-М7(4); 

!define the correlation between residuals (定义 因子 之 间 协 方差 /相关 ); 

R1 WITH МІ; R2 WITH М2; R4 WITH M4; 

R5 WITH М5; R6 WITH M6; R7 WITH M7; 

OUTPUT: STDYX;tech4;tech1; 
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